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Size-dependent effects in solutions of small metal nanoparticles 
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Abstract. A new theoretical approach for the calculation of optical properties of complex solutions is 
proposed. It is based on a dielectric matrix with included small metallic inclusions (less than 3 nm) 
of spherical shape. We take into account the mutual interactions between the inclusions and the quantum 
finite-size effects. On the basis of the effective medium model, TDLDA and Kohn-Sham theories, some 
analytical expressions for the effective dielectric permittivity of the solution are obtained. 

PACS. PACS-key 78.67.Bf - PACS-key 73.20.Mf 

1 Introduction per, we show influence of the quantum-size effects on the 

optical properties of the solutions (matrix disperse sys- 



The optical properties of the solutions with small metal 
particles have attracted much attention in recent years. A 
number of theories have been proposed for describing the 
anomalous far-infrared absorption enhancement in this 
kind of systems 1 . Most authors agree that a direct ac- 
count of the interaction between particles in such systems 
has a significant influence on the behavior of the absorp- 
tion spectra of these systems 2 . In many previous works 
mutual interactions between particles were taken into ac- 
count [2]- [5], though until now the influence of quantum- 
size effects for interacting small metal nanoparticlers on 
the absorption spectra remains uninvestigated. In this pa- 



tems) based on dielectric matrix and embedded in small 
metallic inclusions (less 3 nm) of spherical shape. We demon- 
strate that level of mutual interactions between metal- 
lic inclusions is increased, especially for the small size 
nanoparticles where quantum-size effects became most pro- 
nounced [9]. The underlying mechanism of enhancement 
is related to the difference in the density profiles of sp- 
band and d-band electrons near the nanoparticle bound- 
ary. Specifically, the localized d-electrons are mainly con- 
fined within nanoparticle classical volume while the wave- 
functions of delocalized sp-electrons extend outside of it. 
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This spillover leads to a larger effective radius for sp- neighboring many-particle interactions, can be expressed 

electrons [10] and, thus, to the existence of a surface layer as 

with diminished li-electron population. [191121 As a result. /i o Ar\ (o) , , /i\ 

^ ' p,(l,2...7V) = p^ ^+2^Py +..., (1) 

in the surface layer, the screening of sp-electrons by d-band 

where is a dipole moment of an isolated particle i induced 



by the external field Eq . The second term in the expansion 
(1) is defined by 



P.,=p(*,j)-p!"\ (2) 



electron background is reduced, leading to the enhance- 
ment of the local field and correspondingly to the mutual 
interactions between neighboring metallic particle. 

The paper is organized as follows. In Section [2] we de- 
velop model of extended Maxwell-Garnett approximation where p{i,j) is the dipole moment induced in particle i in 
[13] which will includes all mutual multipolar interactions the presence of a second particle j . The second term in 
between neig hbor metal inclusions [5 . In Section[J]we cal- expansion (1-2) defines pair interactions between particles 
culate absorption spectra for the solution of Au nanopar- in the given mixture, the third three- particle interactions, 
tides. Discussion of our numerical results is presented in and so on. We have restricted ourselves to considering the 
the same Section. In Appendix we propose original way for case of pair interactions between particles, because the 
derivation of quantum multipolar polarizability of metal inclusion of the third term will render the system mathe- 
nanoparticle. Section H concludes the paper. matically hard to treat. As it is follow from the Maxwell 

equations, in the electrostatic approximation the average 

2 Extended Maxwell-Garnett approximation Lorentz field is connected with average polarization in the 

system through the Green function equation 

Consider a system of N non-overlapping spherical metallic 
particles randomly distributed in the background media 

and interacting with each other and with an incident plane From the other side, the averaged macroscopic polariza- 



(F(r)) = Eo + I G(r - r')(P(r'))dr' (3) 



monochromatic wave, whose wavelength A is assumed to tion is linked with averaged Lorentz field in the solution 
be much larger than the average size of the particles in the as 

cluster (i.e., electrostatic approximation) [T3]. The parti- ^p^ — 3 e — Em 

An e + 2sm 



cles have approximately the same radius R and embedded 
in a dielectric continuous matrix with dielectric constant 
Em ■ Due to the influence of the external field as well as 
the interaction between particles in the solution, all parti- 
cles in the system will obtain some polarization. The total 



where e is a effective dielectric permittivity of the solution. 
Since we have restricted ourselves with first two terms of 
expansion (1), the averaged macroscopic polarization of 
the solution becomes 



^ . , , . , (P(l,2...iV))=pf) + l^p,,,(r„r,)<P(r.,r,)dr, (5) 

ettective dipole moment of particle accounting mutual ^ 
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Fig. 1. Two metalhc nanoparticles in the external field Eo, 

where y is a total system volume and <P{ri,rj) is a pair 
correlation function of the particle distribution in the so- 
lution. 

In order to find expression for the pair dipole polariza- 
tion lets us solve two particle problem in the electro- 
static approximation. 

First, we denote by r.^ the center of the sphere i and by 
r an arbitrary point in the medium. The potential inside 
sphere i can be expressed in the form: 

^« = -Eo J2 - ^^I'^irnir - n), (6) 

Ira 

where r — ri is the unit vector along the direction r — ri 
and Yim{r — r.i) are the scalar spherical harmonics (SSH) 
so I is the polar order. The potential outside sphere i can 
be expressed in the form: 



Im 

EoY.Bl:l\r-n\-^-'Y,rn(r-ri) 



(7) 



where the first term is the potential of the external field, 
the second term is the potential created by the particle i 
at the point r , and the third term is the potential created 
by the rest of the particles of the solution at this point. 
Breaking up the potential inside the sphere into compo- 
nents in the two regions and applying standard boundary 
conditions at the three boundaries, plus use of 

dim = -Slid — {V2(5mOU2 + iiSml + <^m-l]Uy 

+ [Sm-l - 5,nl]Ux} (8) 

where u = Eo/i?0: leads to the coupled set of equations 

oo 

^T/;r'Xz-„ = ,5/i; 1^1,2..- m^-l,...l, (9) 



i'=i 



where 



T, 



(m) 



U21+1 { l + V 

-Sw - {~ir{R/n,y+' +1 I (10) 



ai 



I + m 

where i? is a radius of nanoparticlc, and Xim = BimR^^''^^'' 
with symbol 

/ 



I 



II 



m\(l — my. ' 



(11) 



Following the classical approach for polarizability of the 
single spherical particle we should apply an expression 
which can be easily derived from the Mie theory [14] ai = 
_R^'+-^Z(e — £m)/{l£ + {I + l)£m)- Although, if we consider 
small size nanoparticles (less 2nm), that might be neces- 
sary to take into account finite-size effects . The role of 
interband screening and sp band spillover on the bound- 
ary of small particle will make a significant influence on 
the general multipolar polarizability. The quantum mul- 
tipolar polarizabily in Eq. (jlOp can be defined in terms of 
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induced density of electrons [9] 6n = Sng + dnd + Srim due 
to influence of interband screening and sp band spillover 
as = + af where 

/■OO 
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l{2l + 
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-,21 + 1, 



dr'r"+Hn^i\r') 



dr'r'^SrSp (r') 



BiiRd,r') 



a\l + \)\mBi{Ry) 



+ 



Z(2/ + l)i?'+iA„ 
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OO 



/o 

J+i 



Bi{R,r') 



la'+^XdBi{Rd,r') 



with coefficients given in Appendix A 

Finally, from the boundary equations plus Eqs. (5-9) 

we obtain the expression for pair dipole moment 

1 

p{i,j) = -Eo J2 - r,) (13) 

m— 1 

We can now derive an expression for p^j as 

Pij = R^[Xion;,u^ + Xii{n.j.Ux + nyUy)]'Eo, (14) 

where n = ry/ry and rij — |ri — rj\. Here coefficients ATio 
and X\\ are obtained from the coupled system of equa- 
tions Eq.(l9]). By selecting position of two particles along 
the axis Z and parallel to external field Eq, we can investi- 
gate dependence of the pair dipole interaction magnitude 
on the distance between two particles (see Fig (2). Starting 
with small size of particle 92e ox R — O.lnm the maximum 
contribution to the absorption spectra will bring the par- 
ticles located on relatively far distances from each other, 
i.e. 10, 20 and 40 nm. That conclusion lets us to assume 
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Fig. 2. Image part of the pair dipole interaction given in 
Ea. (|14|) in Z direction vs distance between two particles. 
The letters "c" and "q" defines application of the classical 
and TDLDA methods correspondingly for particles with 92e, 
(12) R = 0.7nm; 556e, R = 1.3nm and 2334e, R = 2nm. 



that further in our calculations we can restrict ourselves 
only with the dipole / = 1 interactions between particles, 
neglecting with higher orders multipoles, which are impor- 
tant for the close distances. After performing integration 
over distance between each pair of particle in Eq. ([5]) with 
pair correlation function ^{r-ij) and taking into account 
expression for the Lorentz field in Eq.(|4]), we get a final 
expression for calculation of the effective dielectric per- 
mittivity of the solution with account of pair multipole 
interaction between metallic nanoparticles [S] - [5]. 
e + 2e,n 1 



e - e„ 
1 



AIR' JO 
where parameters 



rlMn,)[f3Hn,) + 2p^{n,)]dn,, (15) 



(16) 



where Ai is given by the I = 1 term in Eq. (jl2p and 
/ = ^R^n is a volume fraction of metallic particles in so- 
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lution with concentration n. The second term on the right 
side of Eq. fTSl) is a contribution of the pair interaction 
between particles. Omitting this term, Ea. (|15p is easily 
reduced to the simple traditional form of the Maxwell- 
Garnett approximation [13] . 

3 Numerical results and discussion 

Below we present results of our calculations for Au nanopar- 
ticles with the same radius ranging from 0.7 to 2 nm in 
a medium with dielectric constant Sm — 1.77. The sur- 
face plasmon (SP) resonance is positioned at 2.35 Ev, far 
away from the interband transitions for gold [18j . Further 
we will consider absorption properties of the given solu- 
tion with metallic inclusions in the SP frequency region. 
The volume fraction of nanoparticle inclusions in the so- 
lution is / = 0.1 and we used the simplest approximation 
for the two particle distribution function 



1 ry > 3i? 



< 3i? 



(17) 



where minimal possible distance between two particles is 
selected to be at least more than one radius or 3i? between 
particle centers. 

To ensure spherical symmetry, only closed-shell magic 
numbers N of electrons 92, 556 and 2334 were used [17 . 
For such sizes, the Au band-structure remains intact. The 
ground-state energy spectrum and wave functions were 
obtained by solving the Kohn-Sham equations for the jel- 
lium model |19| with the Gunnarsson-Lundqvist exchange- 
correlation potential ^D]. These results were used as in- 
put in the numerical solution of TDLDA system A(13). 
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Fig. 3. Calculated real part of the medium effective dielec- 
tric constant in dipole approximation for solution with Au 
nanoparticles with 92, 556 and 2334 electrons at volume frac- 
tion / = 0.1. The solid line corresponds to the calculations 
based on TDLDA method, and the dotted corresponds to clas- 
sical approach. 




Fig. 4. Calculated image part of the medium effective di- 
electric constant in dipole approximation for solution with Au 
nanoparticles with 92, 556 and 2334 electrons at volume frac- 
tion / = 0.1. The solid line corresponds to the calculations 
based on TDLDA method, and the dotted corresponds to clas- 
sical approach. 

The effective radius of nanoparticles was taken in form 
R = r^N^^^, where parameter r^, defined as radius in 
atomic units, is equal to 3 a.u. for Au. On Fig [3] and 
Fig m are shown calculated real and image parts of ef- 
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fective dielectric permittivity of given solution obtained between nanoparticles in solution and on the optical prop- 

from Eq. (jl5[) . Here we can compare absorption [Im{i)) erties of this solution, 
computed with use of TDLDA quantum polarizability a'l 

(layer with sp-free electrons is included with parameter Derivation of quantum multipolar 



polarizability 



A = I, see derivations in[9j) with absorption computed 
classically with polarizability ai. For all three cases (a)- 

(c) on Fig Eland Fig HI results based on TDLDA display The fuU self-consistent potential <f (r) = V'(r) + Vo(r) sat- 

higher level of absorption in the system, especially for the isfies Poisson equation 
solution with small particles (like 92 e). In such case of 
small particle size the difference between two magnitudes 

«'(r)^^o(r) + e2|dV^^, (18) 

can reach almost 2 times and significantly reduces with where ^o{r) = -eEor is potential of the external electric 
increase of particle radius. This fact allows us to predict field, the induced density is comprised of sp-hand, d-band 
that in the solutions with smaU particles the level of ab- and medium contributions, SN{r) = SNs{r) + SNd{r) + 
sorption can be higher than it was expected before based (5iV,„(r). The induced density of sp-band electrons is deter- 
on the classical calculations. mined from TDLDA equation presented below and contri- 
butions from d-electrons 6Nd{r) and medium 6Nm{r) are 
4 Conclusions back to 'P{r) as 

We proposed a theoretical method for calculation of opti- SNd{r) ~ V [xd(r)V!f'(r)] , 

cal properties of the complex solutions based on dielectric SNm{j^) = V [xm(^)V!Z'(r)] (19) 



where Xd{r) = 0{Rd — r) is the interband susceptibil- 
ity with the step function enforcing the boundary condi- 
tions and, correspondingly, Xmir) = 9{r — R) is the 
susceptibility of surrounding medium. The induced den- 
sity of sp-band electrons is given as 



matrix e,„ with included in small metallic inclusions of 
spherical shape with account of the dipole interaction be- 
tween inclusions and quantum finite-size effects. We shown 
that optical properties of given solutions are very much 
dependent on the proper account of microscopic effects, 
in particular, electron screening in the surface layer of 

47r(57ig (^) 

nanoparticles. The exit of delocalized s- electrons out of 5Ns{y) — ^ ^'y2i+i ' 

the classical boundary leads to the additional local field . ^- ^ ^ ^i-oi ^ ^ _li r 

integratmg by parts, Lvq. ([18]) takes the form 

enhancement and as a consequence to enhancement of r xat / i\ 

the surface plasmon of nanoparticle. That correspondingly J l'"^ ~ ^ \ 

makes influence on the character of the dipole interaction 47r J |r — r'| I \ I 



An 
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l{ed - l)a'+V;(^d) + {le,n + + l)MR) 

= i2l + l)MR) (27) 



e„, 1 / ,3 /„/ 1 



dVV— — — • V'6l(r-i?)!f(r'), (21) 



where e(r) = (e^, 1, Cm) for r in the intervals [{0,R), 
{Rd,R), {R,oo)], respectively. We expand and SNg in where a — Rd/R- Substituting ^((i?^) and V'i(-R) back into 
terms of spherical harmonics and obtain, 



Eq. ((241) . arrive at 



e{r)i}i{r) = '0o(r) 
/ drV'2B;(r,r')5A^s(r') 
i?^9fl.,BKr,i?<j)^,(i?d) 



where 



Bi{r.r') = 



+ - 



47r 



47r 



47r 



R^dnBiir, R)MR), 



(22) 



(23) 



21 + 1 

is the multipole term of the radial component of the 
Coulomb potential. The above equation can be simplified 
to 



21 + 1 

Cm - 1 



l3i{T/Rd)MRd) 
Pi{r/R)MR), 



21 + 1 

where we introduced a shorthand notation 



(24) 



Mr) = Mr) + (? / dr'r'^Bi(ry)bN,(r'), (25) 



and /3/(r/i?) = H^R^dRBiiR.r) is given by 



e{r)'>pi{r) = -tpi^r) 



-Pi{r/Rd) — MRd) -{1 + i)a'KnMR) 

T] L J 

+Pi{rlR) — \MR) - la'+^XdMRd)] , (28) 
7/ L J 

where the coefficients A are given by 

(d-l 



a = Rd/R, Ad 



A. 



l€d + l + l' 

T^^l^l{l + l)a''+^\dK 



(29) 



Separating out (^n^-dependent contribution, we arrive 
at TDLDA equation Ref. 

5n^p (r) = / dr'r'^ni^^ (r, r') [^'^ (r') + Sw^''' {r') 



+ / dr'r'^nf\ry) 



dr"r"^A{r'y')5n'i\r") 



+Fi(r')<5n«(r') 



(30) 



where 7Ts(r,r') is the polarization operator for noninter- 
acting sp-electrons, t4'[f^(?'')] the (functional) deriva- 
tive of the exchange-correlation potential and n{r) is the 
ground-state electron density. The latter is obtained in a 
standard way by solving Kohn-Sham equations. Here we 



A(x) = e{x -1)-{1 + l)x^e{l - x). (26) define w'i\r) = r^/e{r) and 



The boundary values of ^ can be obtained by matching 

V'(r) at r = Rd, R, 

{ltd + I + l)MRd) + {1 + l)fl'(em - i)MR) 
= {2l + l)MRd) 



Sw^'\r) = [-Pi{r/Rd)Xda\l - (I + l)X„^)/v 

+ Pi{r/R)X,n{l-la^'+^Xd)/v 
Swi^^r) = [ dr'r'^A{r,r')dni^\r') (31) 



Please give a shorter version with: \authorrunning and \titleruiming prior to \maketitle 



where kernel A(r,r') is a renormalized Coulomb potential 
^{r) due to d-band and medium expressed as 



e(r) 



Bi{ry)-I3i{r/Ri)^ 



-a\l + l)\mBi{Ry) 



+ Pi{r/R) 



BiiRd,r') 
Bi{R,r') 



-la'+^\iBi{Ra,r') 



.(32) 



The expression for quantum nanoparticle polarizabil- 
ity a'^ = ad + cts, can be obtained from the asymptotic for 
r > i? 



{21 + l)e„r'+i 

An 



{21 + l)e„r'+i 



Old, 



as 



(33) 



with 

ad = 



Xda^'+\1 - X^) - X 

pOO 

Jo 

Jo 



+ 



AlTT] 



l{2l + l)i?'+iA„ 



Bi{Rdy) 
-a'{l + l)XmBi{R,r') 



AlTT] 



/ dr'r'2fa«(r' 
Jo 



-la'+^XdBi{Rd,r') 



Bi{R,r') 

(34) 
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